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The liquid-gas phase transition is studied in a multi-component nuclear system using a local 
Skyrme interaction with Coulomb and surface effects. Some features are qualitatively the same as 
the results of Miiller and Serot which uses relativistic mean field without Coulomb and surface effects. 
Surface tension brings the coexistance binodal surface to lower pressure. The Coulomb interaction 
makes the binodal surface smaller and cause another pair of binodal points at low pressure and large 
proton fraction with less protons in liquid phase and more protons in gas phase. 
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O ■ I. INTRODUCTION 

■ 

The liquid-gas phase transition in nuclei was first studied using a Skyrme interaction and focussed mostly on one 
component systems of just nucleons even though expressions were developed for two component systems of protons 

1 ^ ■ and neutrons jlj. The phase transition aspects are considerably easier to study in one component systems rather 

than two, where for example one has to deal with separate proton and neutron chemical potentials for charge and 
nucleon number conservation. Initial fragmentation models also mostly dealt with one component systems. Such 
fragmentation studies gave the first evidence for the liquid-gas phase transition in nuclear systems ^|J3|. Since then, 
the liquid-gas phase transition has been extensively studied experimentally and theoretically. Several reviews exist 
J> ' on this topic |^-^|. 

Because of the two component nature of real nuclear systems, an analysis of liquid-gas phase transition in these 
systems is important. Some preliminary results were reported in Ref. and a very detailed study was done by Miiller 
and Serot Q who used a relativistic mean field model to develop the main thermodynamic properties of asymmetric 
nuclear matter. One interesting new aspect of two component system compared to one component system is that the 
phase transition is a second order transition in their approach. The importance of the number of components on the 
order of the transition was pointed out by Glendenning || . Another interesting aspect of two component systems is 
""^ ' the possibility of having different proton-neutron ratios in the liquid and gas phases because of the symmetry energy. 
The study of nuclear systems with arbitrary proton-neutron ratios is important for radioactive beam experiments and 
1 1 ■ in astrophysical situations as in neutron stars. Because of the extra degree of freedom associated with varying proton 
fraction y in the two phases, the phase diagram has a higher dimensionality and now becomes a surface in pressure 
P, temperature T, and proton fraction y or nucleon density p. For one component systems the phase diagram is 
represented as a binodal curve of P versus density p or volume V, whose end points at fixed T give the liquid and gas 
densities. The Maxwell pressure versus T in one component systems is a line that terminate at the critical temperature 
T c . For two component systems, the binodal surface associated with phase coexistance in (P, T, y) now contains a line 
of critical points for different values of y, a line for the Maxwell pressure versus T for a symmetric system at y = 0.5 
which has a the smallest P and is the same as for a one component system, and a line of maximal asymmetry. At 
fixed T, P versus y of binodal points form loops and (dy / 'cLP)t = gives a point of maximal asymmetry or smallest 
proton ratio on the binodal surface. 

In this paper we extend the initial study of Ref. [0 using Skyrme interaction in a similar way as done in Ref. . 



Some features are qualitatively the same as in Ref. |7fl, but quantitatively differ because of the different interaction. 
Our equation of state based on the Skyrme interaction with Coulomb and surface effects has some feature also not 
present in Ref. Q which will also be discussed. According to the results of Ref. ]?J which has no Coulomb interaction, 
the liquid phase has higher proton fraction than the gas phase of mixed phases in the coexistent state. A real nucleus 
has less protons than neutrons due to Coulomb interaction. For a given A = N + Z stability is determined by 
Coulomb and symmetry energy effects. Since a stable finite nucleus has zero internal pressure while the gas phase 
having positive pressure we also need to consider surface effects. In this paper we consider the effects of Coulomb 
interaction and surface tension by considering uniform spherical finite nuclear system. 

In Section II, the main equations for the thermodynamic properties of hot nuclear matter in mean field theory are 
developed as a function of density p, temperature T, and proton fraction y or neutron fraction (1 — y). These include 
the pressure and chemical potentials, both neutron p n and proton p p where phase equilibrium requires equality of 
the proton chemical potential between the two phases and similarly equality of the neutron chemical potentials at a 
given temperature and P. At fixed T and P, p and y are not independent and we use this connection to simplify 
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the analysis. Section III contains the results of calculations performed using the equations developed in Section II. 
Conclusions are given in Section IV. 



II. PHASE TRANSITION IN MEAN FIELD THEORY 



For phase transition we look at the pressure P and chemical potentials p q for each component of multi-component 
system (neutron and proton for example) as functions of temperautre T and the densities p q of the constituents. 
These quantities can be obtained once we know the total energy functional £ as a function of the densities p q at a 
given temperature T. 

At a given temperature T = 1//3, the constituents are distributed in phase space according to the Wigner function 
/ as 

f(f,0i=^2f 9 (f,0), f g {r,p) = ^W,P) = ^ ^-L) + 1 (1) 

9 

where the spin degneracy 7 = 2 and e q and p q are the single particle energy and the chemical potential of particle 
type q. Then the particle density p becomes 



P(r) = E^( f )' MO = / d 3 pf q (r,p), (2) 

A = J2 N « = J <?rp(f), N q = J <?rp q {r) = J d 3 r J d 3 pf q (f,p) (3) 
and the total energy is 

E = J d 3 r£(r) = J d 3 r J d 3 p^-f(f : p) + J d 3 rU(f) = J d 3 r [£ K {f) + U(r)] (4) 

with the potential energy density U(f) and the kinetic energy density £ft-(r). These, in turn, give self-consistent 
equations for p q (or pFq) in a mean field theory for fixed T and N q , 

9 Sf 9 6f q (r,p) ~2m + 6f q ~ 2m + U ^ (5) 

2 

N = € q\p=p Fq = u i(r,PFq) (6) 



where we define pF q in Eq.(Q) and consider it as an effective Fermi momentum at T of particle q, which is the 

SU 
Sf q 



momentum associated with the chemical potential, and u q = W- is the single particle potential of particle q which 



may be momentum dependent in general. 

The pressure P can be defined dynamically from the total momentum conservation law, 4r [J d 3 r J d 3 p pf] = 

— J d 3 rV r - li— 0, using the Vlasov equation 

^ + (V p e q ) ■ (V r / 9 ) - (V r e g ) • (%f q ) = (7) 
or more generally from hydrodynamic consideration of TDHF in phase space [ pT| 

vv n = -± [/ dVE«^l] = -E/ ^(f) 

= E / d3 PP ' [(V P e g )/ 9 ] + E / S PP ■ (V r e q )f q (8) 
q J q J 

where p is the unit vector in the direction of p. Using (V r e g )/ g = V r (e 9 / 9 ) — e g V r / 9 = V r (e 9 / g ) — V r £, the dynamical 
pressure tensor Uy becomes 
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n y = E / d3 PP* ( V P e «) /« + ^ / d3 pY, e vfi - £ 

1 L q . 

= E / d3 pp* v 



«7 



f q -u 



For a momentum independent potential, this becomes 

SU' 



E 



9 <5p g 



The diagonal elements are 

p = n« = E / + E |^ - u = Pk + E u iPi - u = Pk + Pp 



(9) 



(10) 



(11) 



where Pr- = J d 3 p^f and Pp = ^2 q u q p q — U are the kinetic pressure and interaction (potential) pressure respectively. 
In equilibrium P can also be obtained by minimizing the total energy as a function of volume holding the number of 
particles fixed: 



p = n„ 



d(E/A) = 2 d(£/p) 
dV 9 dp 

for a single component uniform system. 

From the distribution f q of Eq. ([!]), the entropy S can be obtained as, 



s = E S( i = J d3rS = J d3r E 5 f 



S a 



1_ 

h 3 



d 3 p 



I d 3 p-p-(V p f q )ln 

j_ 

h 3 



/,ln/ g + (l-/,)ln(l-/ 9 ) 



1-A 



rfVgP- {V p f q )j3{n q -e q ) 



using V p • (pg) = 3g +p ■ V p g. The iS 9 becomes, by partial integration |l2| ], 



(12) 



(13) 



(14) 



S q = (3 J d p-p- Vp [(/iq - e 9 )/ g ] ~ P j d P3P • v p(^ 9 - £<?)/? 

= /? / d 3 P e q f q + (3 J d'p^^f, - (3p q J d 3 P f q (15) 

In equilibrium thermodynamics, the thermodynamic grand potential f2, the Hclmholtz free energy P, and the Gibb's 
free energy G are 



J d 3 rP= ~ In Tre-K*-^* 1 *^ = F-G = E-TS-J2 



n = 


-j d 3 rP 


F = 


J d 3 rT = 


G = 


/AD- 



PqNq, 



•g = ImTre^^ = 5>,tf 9 = /^E 



(16) 

(17) 
(18) 



3 



Comparing Eq.(p"5|) with Eq.(|^) for diagonal elements using p^VJ, = • V p (isotropic condition in momentum, 
TS = £ + P-Y, W q =£-T = £k + Pk - ~ u i)Pn ( 19 ) 

f = J2 npi - p = Yfa - u ^ p " - Pk + u > (20) 

1 1 

The entropy can also be found from dQ = TdS = CdT; 

The pressure and the chemical potential are also related to the free energy F as 

P— — aF \ — — I /oca 
av \t,a ~ ov It, a ' v °> 

OF I 

The specific heat capacity is given by from the entropy per particle, S/p, 



P=m\rv (24) 



(25) 



To study the caloric curve or the specific heat we look at the enegy per particle £ / p and the entropy per particle S/p. 

For multi-phase multi-component system, the phase transition of each component can occur at different conditions 
such as temperature or pressure. However in general we can treat all the different possible combinations of phases 
of each component in the multi-component system as different phases of the system. (As an example suppose we 
have a system of particles of type p and particles of type n, each having a liquid-gas phase transition but at different 
temperatures. Then the system can be in one of the following four phases; a phase of liquid p and liquid n, a phase 
of liquid p and gas n, a phase of gas p and liquid n, or a phase of gas p and gas n.) Then, to separate each phase of 
the multi-component system, we can use the volume fraction A, of i-phase of total volume V which depends on T, 

Xi = Vi/V with = 1, (26) 

i 

i 

0(p q ,T)=J2WT,pi) (28) 

i 

where O is any obserbable per unit volume. Within the spinodal instability region, there is no equilibrated phase. 
Two or multi phases can coexist when the pressure P and each chemical potential p q are all the same among these 
phases, i.e., 

Pi = Pj, (29) 
H* q - pP q (30) 

with different values of p q and p q for all i and j. At the critical point 

dP d 2 P , s 

dp q dp 2 



1 

dp, q _ d 2 p q 



dp q dp\ 

Spinodal instability occurs when is negative. 



= (32) 



4 



Once the potential energy U in Eq.(Q) is known, then the possibility of a phase transition of the system can be 
studied using Eqs.(|]) - (p5|). The potential energy U determines e q and p q and the potential energy part of E and P. 
Then for fixed T and N q , the Wigner function / and pp q are determined and thus the kinetic terms of E, fi q , and P. 
Using these results, the entropy S and cp can be determined. Relativistic mean field theory is used for the interaction 
in Ref. [7]. The role of the Coulomb interaction and surface energy or the finite size effect of a nucleus are neglected 
in Ref. 7] but will be included in the approach developed here. Their result shows that the neutron evaporates first 
as the energy of the system increases leaving more charge concentrated in liquid phase. In this approach the coupled 
equations of nucleon and mesons lead to a highly nonlinear system and thus the equation of state can be obtained 
only through self-consistent iterative way. These problems may be simplified by using a nonrelativistic zero-range 
Skyrme interaction of 

V12 = t (l + xoP'Wi - r 2 ) + |(1 + x 3 P°)p a (^^p-)6(r 1 - f 2 ) (33) 
For a nuclear system of proton (p p ) and neutron (p n ), this gives the local potential energy density as 

u( Pq ) = |(i + |V - |4 + *») E p\ + + f )p a+2 - |4 + *»K E p\ + c p pI + c ^ ( 34 ) 

Here Cp 13 = ^-e 2 R 2 and C s p n = 47rfi ^ g "(p) — ^yl'/'^ p 2 ^ 3 when we approximate the Coulomb and surface effects as 

coming from a finite uniform sphere of radius R — tqA 1 / 3 with total charge Z (Uc — I e R ^ / ) ■ The typical values for 
the force parameters are given in Table ||. For a symmetric nucleus, N = Z , p q = p/2, and thus 

U( P ) = \t QP 2 + ^hp a + 2 + CpPp 2 p + C sP v (35) 

This potential enegy determines the interaction dependent terms of £, P, e q , and p q which depend on densities without 
explicit T dependence. 

For a momentum independent potential energy as in Eq.(|34|), e q — fi q = (p 2 — p 2 Fq )/(2m) is independent of the 
potential and 

mp) = ^mf), f^p) = e 0^- P % q )/(2m ) + 1 ( 36 ) 

Thus we can evaluate the kinetic terms in £ , P, and p q which are functions of T and pp q . Defining Fermi integral 
F a (v) 



2^_\ a+1 r 00 2 P 2a+1 d P 

l + e *-V< V47tF/ Jo 1 + e 0P 2 /2m-r, q 



FM=I T-^d-=[T^) I , Z^Z^ (37) 



Vg = P (fig - U q) = l3PF q /( 2m ) = PFq/ ( 2 «T) = In Zg, (38) 



A = y 2-kTi 2 /mT (39) 
we can write, for f(r,p) — f(v,p), 

J *pW>& = j?J ^ v^J/^+i = X ~ 3 ^ M > (40) 



Pq 

' ^ ,4: ' ! 



6n 2 ^ 2/3 



^q^\p Kq = jd 3 py q{ r,p ) ^2,JSp 



3 



2m e P(p 2 -p% q )/(2™) 



m ^F 3/2 ( Vq ) = --J=\-*F 3/2 ( Vq ) (42) 

Here ep q is the chemical potential at absolute zero or Fermi energy and pF q is the effective Fermi momentum at T. 
The particle number N q = J d 3 rp(r) determines the effective Fermi momentum pp q (r) or r\ q at T, in terms of density 
Pq(r), 
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r, q (p q ,T) = P( Pq - u q ) = /jg = ^(||a 3 p 9 ) 

For multi(two)-component systems with potential energy given by Eq.(j3j), for a given p q (or pjp q ) and T, 
Pq(Pq,T) = Tij q (p q ,T) 

+to(l + f )p + |(1 + f )(a + 2)p^ - |(1 + x 3 )ap a+1 

-t„(~ + x Q ) Pq + + x 3 )(a - l)2p> - ^|(~ + x^ap^pl 



+C/3p f3 - 1 p 2 p + 2Cp p Pp 6 q , p + v C sP 



P(p q ,T)='£~£ Kq (p q ,T) 



(43) 



(44) 



-|(~ + «o) XX - ^(5 + *»)(« + 



+C(/3 + l)Ap + C s (7?-lK, 



(45) 



|(1 + ?)P 2 - |d + -0) + |(1 + f)p a+2 - %{\ + x 3 )p a ^2p; 



+CpPpl + C sP \ 



TS(p qi T) = X »£Kq(Pq,T) - X]0"9 ~ ?' 



(46) 

(47) 



dQM T ( dS /p\ 

cp(p g ,T) = — — =T — — 



dT 



\ dT 



(48) 



Once we evaluate Fi/ 2 {ti) an d ^3/2 C^?); or more directly 77 = ^^(x) an d ^3/2 C 7 ?)) we can evaluate various thermody- 
namic quantities in terms of and T. 

For low temperature and high density limit, A 3 p large, i.e., when the average de Broglie thermal wavelength A is 
larger than the average interparticle separation p -1 / 3 , we can use nearly degenerate (Fermi gas) approximations E3 ] 
for F1/2 to obtain 



12 V^Fq 



£K q (Pq,T) 



2m 
27 



L , , _ 7T 2 TO 2 / 7 \ 4 / 3 y2 2/3 



A, 



3h 4 \6tt 2 J 



Pq' 



(49) 



^ 3 F 3/2 (Vq) 



■ Pq £ Fq 



5tt 2 



T 



12 \e Fq 



3hr /6tt 



IOto \ 7 



2\ V3 



57T 2 m 2 / T \ 4 /3 



3^ 4 V67r 2 y 



(50) 



In the other limit where A 3 p is small, we have a nearly non-degenerate Fermi gas (classical ideal gas) and the resulting 
equations are given by an ideal gas in leading order with higher order corrections E3] as 



Vqiftqt T ) = P(P-q ~ U q) = m 



2 2 

£Kq(P q ,T) = -P Kq = T^PqT 



7 V 2^: 



P 9 A 3 



2V2 7 

1 PqA 3 /l 2 



2 5/2 7 



a In 
f^ 3 



V 7 
2 



1 

2V2\ 7 



35/2 J V 7 



(51) 
(52) 



G 



For a nuclear system with protons and neutrons with the interaction given by Eq.(Q), the non-degenerate Fermi 
gas limit of Eqs . (]o~l|) and (|5^ ) leads to the following set of equations 



Pq{p,V, T ) = Tln 



A 3 

- ) ^ 



T 



2^2 V 7 



~ Pi 



+Ml + f)p+!(l + f)(a + 2K+ 1 



M\ + *v)p q + %{\ + x a ){a - l)2p a Pq - + x 3 )2ap a - 1 p\ 



+C(3p p - 1 p 2 p + 2CpPp p 8 q , p + vCspi- 1 



Tin 



7 

7 



-)(^ i' 2 n- 11 (It 



'2V2 V 



(a + 2) 



£ 
2 

a+l 



£±(2y-l)£ +7*oP=F( J+a:o)to(2l/-l)(S 



a(2y-l) 2 f£) ±(2y-l)(£)p 



+ 7 C[/3 + 2(l±l)]/+ 1 +C* 



(/3+l±l)(2y-l)(£) /9 + /5(2y-l) 2 (£ 



7,-1 



+rjC s p 

P(p,V,T)=Tp + ^ r2 ( X3 



7 



*3 

12' 



£3. 
2 



+ -0) E ^ - 75(5 + x ^ a + E ^ 



(53) 



+C(/3 + lW + C.fa -!)/>•> 



2 ' (a+1 v+ 2 



T 



16 



a + l 
6 



2V2 V 7 
1 



G) (f ) 



*3 7T + *3 p" ~ 



\(P + V C pe +2 + (r ) -l)C s p r > 
kT /A 3 



2^2 V 7 



-\-(j3 + l)CpP 



tey-i f f? 



+ (/3+l)Cp^ 1 (2j/-l) (0, 



(54) 



+CA 2 + C s p" 

3^ 3 2 1 , o 3 T /A 3 
2^+8^ +16^ + + 2i7f(T 



+Cp^\2y -!)(£ 



x 



1 



3 feT / A 



TS(p,y,T) = -7>- T£p g ln — 



2 



A 3 



■Pa 



— — y In 

2 y 



T 



2V2 V 7 



A 3 



/A^ /EgP 



-(l-y)ln — (l-y)p 



(2y-l) 2 (J 



(55) 



7 



T 



A 3 \ [l + (2y-l) 2 



(§)" 



2V2 V 7 J 2 
Here, for the proton density (p p ) and neutron density (p n ), we defined, 

P = P P + Pn, P3 = pp - Pn = (2y - l)p, y = p p /p, 

Pp = Hp + Pa) = yp> Pn = Up ~ ^3) = ( : - 2/V, 

E g ^ = Kp 2 + Pa) = fl+(2 r 1)2 V = [1 + 22/(2/ - 1)]P 2 , 
E 9 P 3 . = ^(P 2 + 3p|) = [1+3(2 r 1)2] P 3 = [1 + 32/(2/ ~ 1)]P 3 



(56) 



(57) 



The ± in p q stands + for q = proton and — for neutron. For £3 7^ —1/2 and 0, the P(p) curve for different vales 
of y at fixed T may cross at some density p which was not seen in Ref. . 

For a constant T and constant P, p and y are not independent. Pressure P of Eg. (|54|) is a second order polynomial 
of (2y — 1) (I), and thus we have, for the range of < y < 1/2, 



1< (22/ -1) 



;i + Xo) * + («±i) (i + .3) hP « - 09 + 1) Cp? - ^ (f 




(/3 + l)<V 



) ^ + m (i + ^ 3 ) ^ -03+1) ^ - ^ (f ) 



^0 



'0 + (¥) S + ^) ^ - (/3 + 1) - ^ (f) 



1/2 



< 



(58) 



for a given density p. (Here + sign is allowed for the case that the first term is negative.) Without the Coulomb 
interaction only the second term in the square root survives and y{p) is a single valued function of p in the range of 
< y < 1/2. The numerator of the second term in the square root is P(p,y = 1/2, T) — P(p,y,T), i.e., the negative 
of the pressure measured with respect to the pressure for a symmetric nuclear system. Notice here that to is negative 
and P(y) > P(y — 1/2) for the potential without the Coulomb interaction (see Eq. (|5~i|) ) . Since we are considering 
only —1 < (2y — 1) < 0, we have conditions of 



< 



< 



Tp 



¥op 2 



(g±i) 

16 



hp 



,a+2 



(/3+1 



•Cp' 3 + 2 + (ry-l)C^ + ^ 5 (^)^-P 



(i + XQ ) tQ + (s±i) (i + 3.3) t3pa _ { p + i) C pP - ^ (f ) 



209 + 1)<V 



; 1 + Xo) to + («±i) (i + X3 ) Hpa _ 09 + 1) - ^ (f ) 



(59) 



to have < y < 1/2 of which the boundaries are related to the pressures at y = 1/2 and y = respectively. For non- 
zero C. there is a back bending in y(p). This behaivor is related to the choice of opposite sign for the square root term 
in Eq.([58|) and double values of y for a given p. To find the point of back bending, we look at the (dp/dy) T P — point. 
Using Eq.(j5"§|) we can find the energy and entropy as a function of pressure P, proton fraction y, and temperature 
T, i.e., £(P,y,T) and S(P,y,T). From these results we may study an isobaric phase transition. Since all the phases 
have the same pressure in the coexistent region, p q (P, y, T) may be used to find the coexistent region. 

At fixed T and P, only one of either p or y is the independent variable. Thus observables such as P, £/p, S/p may 

have a discontinuity in T or y when \§xj or (fa) diverges. We can study the behavior of thermodynamic 
quantities at a fixed P using dP = from Eq.(p4|), 



dP = 




8 



*0(l + ^V + |(l + ^Hn+ l)(n + 2. / /^- > 



X3> 

2 



-*o(~ + *o) (E ^) - §(5 + ^)(« + !)(«+ 2)p Q ^ ^ 



-COS + 1)09 + 2)fPf$ + C s (v- 1)VP V ] } 



to ( 2 + 



v 

a + 1 
6 



dp 
P 



-(/3 + l)C7p' 3 (0 2 |4d 2/ 



1 1 



A 3 



1 1 



fh-" 2 (f)1h 



2 2\/2 V 7 / V 2/ 2 2V2 
T + 2b 2 p +(a + 2)b 3 p a+1 + ^= (£\ (0 + 09 + 2)& c / +1 + ^ - lJC.p"" 1 



*0 ( 2 + ^0 



x (2» If (f 



'11 ( 2" + Xo 



a + 2 



7 



(/3 + 2)(/?+l)C/(2y-l)(0}dp 



-(/9+l)Cp 



P(P 



a + 1 



4dy 



-.(l + ^^-^ + DC/-^^ 



(%-D(f) : 



Thus we have, for a fixed P and y, 

T /A 3 



<9p 
ST 



7>- 



2V2 V 7 
T /A 3 



2V2 V 7 



toil + yV + g(l + f )(a + l)(a + 2)p«+ 2 
-*<4 + *o) (E P^ - ^(5 + x s )(a + l)(a + 2)p« ^ p 2 



_P 

T 



+C09+l)09 + 2)/pg + a.(» 7 -l)» ? p'']} 1 
1 1 /A 3 \ /p\2 1 l ( A 3 



(I)' 



2 2 V / 2V7/V2/ 2 2V2V7 
T + 26 2P + (a + 2%p a+l + ^= {^-\ (|) + 08 + 2)6 C / +1 + 77(77 - lJCp"" 1 



•i'n 



a + 2 



/? + 2 



09 + l)C/ 



2\/2 V 7 



A 3 



x (2y - I) 2 (£) + 09 + 2)03 + l)C/(2y - 1) (f ) } * 
Using this we get the specific heat capacity, from Eq.(E6 



rp ( dS / p\ 



T 



5 3 , /A 3 
2 + 2— yXn \- yP 



/A 3 



(l-y)lnl — (l-y)p 



1 T /A 3 \ [l + (2y-l) 2 ] /p 



2 2^2 V 7 



(!) 



(60) 



(61) 







-T 



T /A 3 \ [l + (2y-l) 2 ] (p 



T 



4 - y In 
T 



2V2 V 7 
'A 3 



r 



<9T 



7 



-VP 



II- ,/)]„( ^-(l-y)p 



1 T /A 3 \ [l + (2y-l) 2 ] /p 



2 2^2 V 7 



(?) 



T 



_/A 3 \ [l + (2y-l) 2 ] ,ps 
2V2V7/ 2 W 

T + 2b 2P +(a + 2)b 3 p a+1 + (^) (|) + (/? + 2)6 C / +1 + 77(77 - l^p"- 1 



to 



Xq 



a + 2\ fa + 1 
IT" 



t 3 



+ Z3 }P 



13 + 2 



(/3 + l)C/ 



T / A 



2V2 V 7 



(2y I) 2 (f ) + 08 + 2)(/3 + l)C/(2!/ - 1) (|) } * 



4Tp-T^p g ln 



A 3 



-Pa 



1 T /A 



2 2^2 V 7 



2V2 V 7 

T /A 3 



'HqP. 



kTp + 



2V2 



7 ?' 



*o(l + y)p 2 -*o( 2 



+|(l + y)(« + l)(- + 2K +2 - 1(^+^3 
+C(/3 + 1)09 + 2) A 2 + C s (ry - i)r,p"] , 



T (ds/ P \ 



r 



A 3 



7 



4- 77 In — yp - (1 - y) In —(1 - y)p 



A 



4Tp-TEp 9 ln 



-P9 



7 



1 T /A 



1 T /A 3 \ [l + (2y-l) 2 ] (P 



2 2^2 V 7 



(5) 



2 2V2 V 7 



This shows that there is a divergence of cp at T with ( 

transition. This happens when 
volume cy has no divergence. 



00, i.e. 



(f) 



(62) 



(63) 

= showing first order phase 



T,P 



which may exist in the coexist region. The specific heat with fixed 



III. PHASE TRANSITION OF A FINITE NUCLEUS WITH COULOMB INTERACTION 

Here we use Skyrme interaction of PRC45 in Table | with j3 = 0, 77 = 2/3, the surface energy Amr^a = 20.0 MeV, 
and R = 6 fm. Fig.|l| shows the pressure P(p) for various y values at T = 10 MeV. With x 3 = —1/2, there is no 
crossing of P(p) curves between different values of y at fixed T which is a behavior similar to that in Ref. 0. However, 
for T = 10 MeV, the lowest pressure occurs not at y — 1/2 but at y = 0.4552 for R — 6 fm due to the Coulomb 
repulsion independent of the surface tension (see Eq.(|54|)). For R = 10 fm, the lowest pressure occurs at 77 = 0.3927. 

For fixed P and T, y and p are not independent. Since Eq.(|54|) is a second order equation for y, y can have two 
values for a given density p in general. Fig.^J shows y(p) for various P values at T = 10 MeV. Here only the < 77 < 0.5 
region is shown. These curves show back bending for large y with |£ = at y — 0.4552 < 1/2 for the parameter 
with Coulomb interaction which is used here. At this y value the pressure is minimum as mentioned above. Without 
Coulomb interaction, there is no back bending in the < y < 0.5 region and ^ = at 77 = 1/2. The discontinuity of 
y(p) at low density for small pressure is expected from the low density region of Fig.|l|. 
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Fig.|3j shows chemical potential p,„(y) for various P values at T = 10 MeV. From this figure we can find binodal 
points by checking the condition Eq. (|32|) . There is a cross between the small density portion and high density portion 
of the curves for low pressure (see the box at the right-bottom corner of Fig.||), which has a discontinuity in y(p) 
shown in Fig.|[ for both p p and p n . The overall behavior is qualitatively similar to the result of Ref. 0. However with 
Coulomb interaction p n < p p at y = 1/2 instead of p n — p p . The existence of these crossings of chemical potential 
curves between high and low density and between proton and neutron allow two pairs of binodal points at low P as 
shown in the box at the right-upper cornner of Fig.^. 

Fig.[| shows the coexistance binodal curve in y-P plane at T = 10 MeV. The figure shows that the major effect of 
the surface tension is a lowering of the pressure for the coexistence curve (compare dash-dotted curve vs dashed curve 
and solid curve vs dash-dot-dot-dotted curve). Inclusion of surface tension may allow for a zero pressure isobaric phase 
transition and may simulate the situation of an equilibrated state of multifragmentation having zero internal pressure 
of stable finite nuclei with non-zero gas pressure as discussed in Ref. Q . The effect of the Coulomb interaction makes 
the coexistence curve smaller. However the more important effect of Coulomb interaction is the appearence of another 
pair of binodal points in the low pressure region (see the expanded figure in the box at the right upper corner). Two 
points of smaller y values have the same chemical potential and two points of larger y values have the same chemical 
potential but which are different (having increased p p and decreased (j, n ) from the other two (see Table [n]). Here the 
points with largest and smallest y value are in gas phase and the points between the other two y values are in the 
liquid phase (see Table |l]) . The binodal pair with higher y values allows the mixture of gas with highest y and liquid 
with the next high y and is a reversed situation from the other binodal pairs and the result of Ref. H . This allows 
for a phase separation into a liquid phase with less proton (lower y value) and a gas phase with more proton (higher 
y value) which is more realistic in heavy ion collision. Two pairs of binodal points meet together with the lowest 
pressure of the binodal curve at one point with y = 0.4552 for R = 6 fm which would be changed to y = 0.3927 for 
R = 10 fm. Without Coulomb interaction, the lowest pressure of the binodal curve occurs at y = 1/2. 

Fig.|| shows T(p) for various P values at y = 0.3 and 0.5. This figure shows that at a low pressure there are two 
points with = Q which make thermodynamic quantities have a first order discontinuity with T variations. 

Fig.| and Fig.@ shows S(T)/p and S{T)/p for various P values at y = 0.3 and 0.5. Comparing with the exact result 
(dotted curve in Fig||), the Fermi gas approximation used here is good for the region which corresponds to £ / p > 8 
MeV for the case without Coulomb and surface effects. The points with diverging slope of £/p and S/p in variations 
with T are related to the zero slope points in Fig j^. If the charge concentration y was kept constant in all phases then 
these two diverging points might be outside of the coexist region and cause a first order phase transition. However in 
the actual phase transition of a heavy ion collision these diverging points belong to the coexistence region and gives 
a second order phase transition as in Ref. M. 



IV. SUMMARY AND CONCLUSION 



In this paper we studied the liquid-gas phase transition in a hadronic system made of protons and neutrons or 
having two components. The approach used here was based on a mean filed theory and we calculated the mean 
potential energy using a Skryme interaction. The work is an extension of a previous study initially begone in Ref. 0] 
for two components and the analysis is similar in some ways to that developed in Ref. J7jj. However, in this paper we 
included surface and Coulomb effects so our results differ quantitavely from those in Ref. . Including Coulomb and 
surface effects introduces some feature also not present in Ref. 

The importance of surface energies in the liquid-gas phase transition was shown in Ref. using an approach 
based on a statistical model of multifragmentation. In the statistical model of multifragmentation, because nuclear 
matter can be broken into a large number of small pieces, the surface energy plays a very dominant role. Here, in a 
mean field theory, the surface energy plays a lesser role since we do not allow for the possibility of multifragmentation. 
Rather, the system is treated as uniform but with varying density. Surface tension here brings the coexistance binodal 
surface to lower pressure which could allow an isobaric transition at zero pressure. 

Inclusion of Coulomb interaction brings a new feature in a mean field approach. The Coulomb interaction makes 
the binodal surface smaller and cause another pair of binodal points at low pressure with a reversed proton fraction 
in the liquid-gas phase, i.e. with less protons in liquid phase amd more protons in gas phase. The value of the proton 
fraction y for the liquid phase of this binodal pair is larger than the y values for the lowest pressure point of the 
binodal curve as can be seen in the box in Fig.[| and we can make this value smaller by using a stronger Coulomb 
interaction. 

In this paper we used a fixed value of R to determine the strength of Coulomb interaction and surface tension. If 
we allow for variations of R with p or T and allow different values for the Coulomb and surface energies, then this 
simple model may be used to study multifragmentation with various size of clusters. 
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TABLE I. Parameter sets for Skyrme interaction pJJl 



Force 


a 


to (MeV-fm 3 ) 


Xo 


t 3 (MeV-fm 3(1+Q °) 


X3 


PRC45 


1 


§Ci = -1089.0 


1/2 


-^>C* 2 = 17480.4 


-1/2 


ZR1 


1 


-1003.9 


0.0, 0.2, 0.5 


13287.2 


1.0 


ZR2 


2/3 


-1192.2 


0.0, 0.2, 0.5 


11041.0 


1.0 


ZR3 


0.1 


-4392.2 


0.0, 0.2, 0.5 


26967.3 


1.0 
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TABLE II. Coexistance points. 



P (MeV/fm 3 ) 


p (fin 3 ) 


y 


(MeV) 


(MeV) 


0.01496 


0.0034324 


0.4526870 


-15.73360 


-13.80920 




0.1138296 


0.4548638 


-15.73360 


-13.80920 




0.1138296 


0.4555570 


-15.59369 


-13.92612 




0.0034324 


0.4577830 


-15.59369 


-13.92612 


0.01498 


0.0034352 


0.4441228 


-15.97430 


-13.61531 




0.1138098 


0.4536912 


-15.97430 


-13.61531 




0.1138096 


0.4567391 


-15.35914 


-14.12939 




0.0034352 


0.4665239 


-15.35914 


-14.12939 


0.01500 


0.0034373 


0.4253351 


-16.50124 


-13.19126 




0.1137651 


0.4511254 


-16.50124 


-13.19126 




0.1137652 


0.4592869 


-14.85386 


-14.56769 




0.0034372 


0.4853333 


-14.85386 


-14.56769 


0.01503 


0.0034421 


0.4226954 


-16.57709 


-13.12985 




0.1137598 


0.4507550 


-16.57709 


-13.12985 




0.1137592 


0.4596993 


-14.77187 


-14.63844 




0.0034420 


0.4883935 


-14.77187 


-14.63844 


\J .KJ lOKJO 




f) 41 41 1 77 








0.1137397 


0.4495787 


-16.81862 


-12.93537 




0.1137389 


0.4608853 


-14.53660 


-14.84241 




0.0034448 


0.4971483 


-14.53660 


-14.84241 


0.01506 


0.0034463 


0.4102457 


-16.92777 


-12.84748 




0.1137306 


0.4490471 


-16.92777 


-12.84748 
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FIG. 1. Pressure P(p) versus p at T = 10 MeV for various y. The dashed curves from top to bottom have y — 0, 0.1, 0.2, 
0.3, 0.4, and 0.5 while the solid curve is for y — 0.4552. 
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FIG. 2. Proton fraction y{p) for P = 0, 0.015, 0.05, 0.1, 0.2, 0.3, 0.4, and 0.5, from top to bottom, at T = 10 MeV. The 
solid curves have Coulomb and surface tension contributions and the dashed curves are without these interactions. 
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FIG. 3. Chemical potential n P (y) and ^ n (y) for P — 0.015, 0.05, 0.1, 0.2, 0.3, 0.4, and 0.5, from top to bottom curve for 
proton (solid curve) and bottom to top for neutron (dashed curve) at T = 10 MeV. Chemical potential increases as pressure 
increase at y = 1/2. fi n increases and fi p decreases as P increases except fi p for P > 0.05 or 0.1 for smaller y. The small box at 
the right-bottom corner is the expanded curve for P = 0.015 MeV with two pairs of binodal points indicated by dotted lines. 
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FIG. 4. Binodal curve at T = 10 MeV. The solid curve is for the case with Coulomb and surface effect, the dash-dotted 
curve is for the case with Coulomb interaction, the dash-dot-dot-dotted curve is for the case with surface effect, and the dashed 
curve is for the case without Coulomb and surface effect. Small boxes at the upper right corner are the expansion of the main 
figure. In small boxes, the region of liquid, gas, and coexistence are indicated by L, G, and C respectively. 
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FIG. 5. Temperature T(p) for P = 0.0, 0.015, 0.1, 0.2, 0.3, 0.4, and 0.5, from bottom to top, at y = 0.3 (dash-dotted curve) 
and 0.5 (solid curve). Dashed curves are for the case without Coulomb and surface effect. 



18 




-20 ' ' 

5 10 15 20 25 

Temperature, T (MeV) 



FIG. 6. Energy per particle £(T)/p for P = 0.0, 0.015, 0.1 0.2, 0.3, 0.4, and 0.5, from top to bottom, at y = 0.3 
(dash-dot-dot-doted) and 0.5 (solid). Dashed curves are for y = 0.5 without Coulomb and surfase effects and dotted line is for 
exact result at P = from Ref . jl3| . This shows that the Fermi gas approximation we have used is not varied for T lower than 
about 5 MeV with high density, but very accurate for higher temperature (dashed curve for P = 0.0 coincide with the dotted 
curve for £/p> 8 MeV). 
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FIG. 7. Same as Fig.^ but for entropy per particle S(T)/p 



20 



